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Abstract: We present a review of the Bielefeld-Dubna activities on 
the multiloop calculations. In the first part a C-program DIANA (Diagram 
t O ■ ANAlyser) for the automation of Feynman diagram evaluations is presented, 

. in the second part various techniques for the evaluation of scalar diagrams 

t^j- ' are described, based on the Taylor expansion method and large mass ex- 

(Nj . pansion. 

y—i • 
>'■ 

C^- , 1. Automation of Feynman diagram evaluation 

cn ■ 

Recent high precision experiments require, on the side of the theory, 
high-precision calculations resulting in the evaluation of higher loop dia- 
, grams in the Standart Model (SM). For specific processes thousands of 

Q^ ■ multiloop Feynman diagrams do contribute, and it turns out to be im- 

possible to perform these calculations by hand. This makes the request for 
p • automation a high-priority task. 

Several different packages have been developed with different areas of ap- 
plicability. For example, FEYNARTS / FEYNCALC Q are MATHEMA- 
TICA packages convenient for various aspects of the calculation of radiative 
corrections in the SM. There are several FORM packages for evaluating mul- 
• r- 1 ! tiloop diagrams, like MINCER ||], and a package Q for the calculation of 

^ ■ 3-loop bubble integrals with a mass. Other automatic packages are GRACE 

[Q] and COMHPEP [||, which partially perform full calculations, from the 
process definition to the cross-section values. 

A somewhat different approach is persued by XLOOPS ||. A graphical 
user interface makes XLOOPS an 'easy-to-handle' program package, but is 
mainly aimed to the evaluation of single diagrams. To deal with thousands of 
diagrams, it is neccessary to use special techniques like databases and special 
controlling programs. In ]?]] for evaluating more than 11000 diagrams the 
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special database-like program MINOS was developed. It calls the relevant 
FORM programs, waits until they finished, picks up their results and repeats 
the process without any human interference. 

All these packages have different efficiency in different domains. It seems 
impossible to develop an universal package, which will be effective for all 
tasks. This point of view motivated us to seek our own way of automatic 
evaluation of Feynman diagrams. 

Our first step is dedicated to the automation of the muons two-loop 
anomalous magnetic moment (AMM) \{g — 2) . For this purpose the pack- 
age TLAMM was developed [||. The algorithm is implemented as a FORM- 
based program package. For generating and automatically evaluating any 
number of two-loop self-energy diagrams, a special C-program has been 
written. This program creates the initial FORM-expression for every di- 
agram generated by QGRAF fj~C)[| , executes the corresponding subroutines 
and sums up the various contributions. In the SM 1832 two-loop diagrams 
contribute in this case. The calculation of the bare diagrams is finished. 
For the purpose of demonstration, we have applied TLAMM to a closed 
subclass of diagrams of the SM which we refer to as "toy" model. Some 
details of the calculation are presented for this case. 

A more general project called DIANA (Diagram ANAlyser) Q for the 
evaluation of Feynman diagrams is being finished by our group at present 
and will also be shortly described below. 



1.1. The toy model 

We considered as a toy example the model involving a light charged 
spinor ty, a photon A^, and a heavy neutral scalar field <£. The scalar has 
triple (g) and quartic (A) self- interactions, and the Yukawa coupling (y) to 
the spinor. The Lagrangian of this model reads (in Euclidean space-time) 

l = ^d^d^ + x -m 2 & - |d> 3 - + i ^ Av _ duAfi) 2 

+— (d^) 2 + # (d + ieA + m) * - y$\j>tf (1) 

where e is the electric charge and a is a gauge fixing parameter. 

There are 40 diagrams contributing to the two-loop AMM of the fermion. 
The calculation of the AMM reduces, after differentiation and contractions 
with projection operators, to diagrams of the propagator type with external 
momentum on the fermion mass shell (for details see ) . After performing 
the Dirac and Lorentz algebra, all diagrams can be reduced to some set of 
master integrals. 
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Fig. 1. The prototypes and their subgraphs contributing to the large mass expan- 
sion in the toy model. Thick, thin and dashed lines correspond to the heavy-mass, 
light-mass, and massless propagators, respectively. Dotted lines indicate the lines 
omitted in the subgraph 7. 



For the evaluation of these integrals we use the asymptotic expansion in 
large masses [12]. For a given scalar graph G the expansion in large mass 
is given by the formula 



F G (q, M, m, e) M ~°° ]T F G/j (q, m, e) o T gW F 7 (<f , M, m\ e), (2) 



where 7's are subgraphs involved in the asymptotic expansion, G/j denotes 
shrinking of 7 to a point; is the Feynman integral corresponding to 7; 
Tq 1)rril is the Taylor operator expanding the integrand in small masses {m^} 
and external momenta {q-y} of the subgraph 7 ; o stands for the convolution 
of the subgraph expansion with the integrand Fqj^. The sum goes over 
all subgraphs 7 which (a) contain all lines with large masses, and (b) are 
one-particle irreducible w.r.t. light lines. 

Individual integrals are specified by the powers of the scalar denomina- 
tors, called indices of the lines. From the point of view of the asymptotic 
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expansion method the topology of the diagram is essential. All diagrams 
of the toy model that contribute to the two-loop AMM can be classified in 
terms of 9 prototypes (we omit the pure QED diagrams). These prototypes 
and their corresponding subgraphs involved in the asymptotic expansion, 
are given in Fig. |l|. Almost all integrals occuring in the asymptotic expan- 
sion of the muon AMM in the SM can be evaluated analytically using the 
package SHELL2 @. 



1.2. The piloting C-program TLA MM 

For the automatic calculation we have created a special piloting program 
written in C. This program 

1. reads QGRAF output; 

2. creates a file containing the complete FORM program for calculating 
each diagram; 

3. executes FORM; 

4. reads FORM output, picks out the result of the calculation, and builds 
the total sum of all diagrams in a single file which can be processed 
by FORM. 

Identifiers for vertices and propagators and the explicit Feynman rules are 
read from separate files and then inserted into the FORM program. Be- 
cause the number of identifiers needed for the calculation of all diagrams at 
once may exceed the FORM capacity, the piloting program retains for each 
diagram only those involved in its calculation. 

All initial settings are defined in a configuration file. The latter contains 
information about the file names, identifiers of topologies, the distribution of 
momenta, and the description of the model in terms of the notation that is 
some extension of QGRAF's. The program carries out the complete syntax 
check of all input files except the QGRAF output. 

There exist several options which allow one to process only the diagrams 

• explicitly listed by number; 

• of a given prototype; 

• of a specified topology. 

There are also some debugging options. 

QGRAF generates the diagrams in symbolic form in terms of vertices. 
Fig. |2| shows one of the diagrams in the toy model with the corresponding 
QGRAF output. 
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*— #[ d2: 
* 

1 

*vx(E2(2) ,e2(-l) ,H(1)) 
*vx(E2(3) ,e2(4) ,A(-3)) 
*vx(E2(-2) ,e2(5) ,H(6)) 
*vx(E2(5) ,e2(3) ,H(1)) 
*vx(E2(4) ,e2(2) ,H(6)) 

* 

*--#] d2: 




Fig. 2. Diagram "number 2". e2 stand for fermion, E2 for anti-fermion, H is a 
scalar, and A the photon. External legs are numbered by negative numbers and vx 
stands for the vertex for the particles in the argument list. 



The Feynman integrand generated by the C-program for this diagram 
reads: 

g VF(mu)= 

E2_H_e2*F_F (-k2+p , q/2 , me2) *E2_H_e2*F_F (kl-k2 , q/2 , me2) * (3) 
A_E_e (mu) *F_F (kl-k2 , -q/2 , me2) *E2_H_e2*F_F (kl , -q/2 , me2) * 
E2_H_e2*H_H(kl-p,0)*H_H(k2,0) ; 

A demo diskette for evaluating the toy model can be provided on request. 



1.3. Project DIANA 

At present we have finished the main part of the more general program 
DIANA (Diagram ANAlyser). For generating Feynman diagrams we use 
again QGRAF. The program DIANA consists of two parts: 

• Analyzer of diagrams. 

• Interpreter of a special text manipulating language (TM); 

The TM language is a very simple TeX-like language for creating source 
code and organizing the interactive dialog. 

The analyzer reads QGRAF output and passes necessary information to 
the interpreter. For each diagram the interpreter performs the TM-program, 
producing input for further evaluation of the diagram like for example (H|) . 
Thus the program: 

Reads QGRAF output and for each diagram it: 
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1. Defines the topology. 

2. Looking for it in the table of all known topologies and distributes 
momenta according to the current topology. 

3. Creates a special internal representation of the diagram corresponding 
to the Feynman integrand. 

4. Invokes the interpreter to execute the TM-program and passes to it 
the necessary data. 

Executing the TM-program provides the possibility to calculate each 
diagram using FORM or another formulae manipulating language, to do 
some numerical calculation by means of FORTRAN, to create a postscript 
file for the picture of the current diagram, etc. 

If we do not yet know all needed topologies, we may use the program to 
determine missing topologies that occur in the process. 

The main goal of the TM-language is the creation of text files. In prin- 
ciple we could have used one of the existing languages, but we want a very 
specific language: it should be powerful enough to create arbitrary program 
texts. On the other hand, it should be very simple and easy-in-use, so that 
even non-programmers can use it. 

Similar to the TeX language, all lines without special escape - characters 
("\") are simply typed to the output file. So, to type "Hello, world!" in the 
file "hello" we may write down the following program: 

\program 
\setout (hello) 
Hello, world! 

Each word whose first character is an escape character will be considered 
as a command. This feature makes this language very easy-to-use. 
At present, we have finished the C-part of this project. 



2. Evaluation of scalar diagrams 

In this second part methods for the evaluation of scalar three point func- 
tions will be discussed. We will concern ourselves with the Taylor expansion 
with respect to an external momentum squared and the large mass expan- 
sion explained above. The efficiency of both approaches will be compared. 



2.1. Expansion of three-point functions in terms of external momenta 

squared 
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Taylor series expansions in terms of one external momentum squared, 



q say, were considered in [14|, Pade approximants were introduced in [15] 



and in Ref. [16] it was demonstrated that this approach can be used to cal 
culate Feynman diagrams on their cut by analytic continuation. The Taylor 
coefficients are expressed in terms of "bubble diagrams" , i.e. diagrams with 
external momenta equal zero, which makes their evaluation relatively easy. 
In the case under consideration we have two independent external momenta 
in d = 4 — 2e dimensions. The general expansion of (any loop) scalar 3-point 
function with its momentum space representation C(pi,p 2 ) can be written 



as 



C(pi,p 2 ) 



oo 

E 

l,m,n=0 



aimn(pi) l {pl) m {PlP2Y 



(4) 



where the coefficients a/ mn are to be determined from the given diagram. 



For many applications it suffices to confine to the case p\ 



n 



0. 



which is e.g. physically realized in the case of the Higgs decay into two 
photons (H — ► 77) with p\ and p 2 the momenta of the photons. In this case 
only the coefficients aoon are needed. 

In the two-loop case we consider the scalar integral (k$ = k\ — k 2 , see 
also Fig. 1) 



C(mi, ■ • ■ ,m e ;pi,p 2 ) 
1 ■/ 



d i k 1 d 4 k 2 



(5) 



WJ* J ((fci+pi) 2 -m^)((fci+p2) 2 -m^)((fe2+Pi) 2 -m^)((fei+P2) 2 -m^)(^-m^)(^-m^)- 

kf in line 4 (with mass 771,4) depends on the topology: for the planar diagram 
we have kf = k 2 while for the non-planar we have kt = k%. 
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Fig. 3. Planar and non- planar scalar vertex diagrams and their kinematics 



With obvious abbreviations for the scalar propagators: Cj the i th scalar 



propagator of (||) with p\ = p 2 = 0, we can quite generally write for the 



th 
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Taylor coefficient [18|: 

(^ 2 ) 2 a 00 n = -— / d^hd^Fn . (6) 

n + 1 J Cl C2 C3 C4 C5 C6 

The numerator can be written as 

n n 

F n = Yl ^ n - v) c r E ^ {n ~ v V' • ^(fci^, • (7) 

^=0 i/=0 

with 

A1AkiM,h)= E EEC'" (8) 

M =0 cr=0 T=0 

where the coefficients h n v ^f T are given by 

= (n + 1) 2 A -™- 1 (n - u)\ (n - v')\v\ v'\- -t- (9) 



2 2 " d 

^E(^-i)ir(n+ 

r i=l 



min(ni ,n— f) min(?ii ,n— u') nj+k 

y y 

^max^n,-,) k=m^n % - V ') P)K*i ~ k ~ p)l(p + J + k - m)\ 

1 



a j \{ Vj -2a)W k \{u' k -2T)\ 



with X = i/ + i/ — 2fi ) p = fj l — <T — T,Vj = u — ni+j,u' k = u' — ni + k,<Tj = 
a — \j — (n — v) + (i — 1)] and = r — [A; — (n — z/) + (i — 1)]. 
For the planar diagram (k t = £2) (@) simplifies and we have 

^(*l,fc) = E C't*!)^'^''*^*')^*! fe)"^- 2 ", (10) 

with 

<C = E E C T . (id 

cr=0 r=0 
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The coefficients a^, are mass independent and have been calculated 
with FORM up to order e 2 (d = 4 — 2e) and stored for the first 30 Tay- 
lor coefficients, i.e. they are given in terms of rational numbers. For the 
non-planar case the situation is more difficult since the storing of the coeffi- 
cients b n J^f T with two more indices is practically impossible for high Taylor 
coefficients. In many cases, however, if the threshold of the diagram under 
consideration is high, only a few (say 8) Taylor coefficients are sufficient for 
a high precision calculation of the diagram under consideration and in such 
a case the direct calculation of the b's in each case causes no problem. 

Finally all remaining integrals can be reduced to bubble-integrals of the 
type 

v ,r f d d k X d d k 2 

Q/Mim|J J {k 2 -m 2 )«(k 2 -m 2 )P{{k x -fc 2 )2_ m 2) 7 ' [lZ> 

or to factorizing one-loop integrals. The genuine two-loop bubble integrals 
are reduced by means of recurrence relations to Vm({m}). This can be 
done numerically for the arbitrary mass case or also analytically for special 
cases like e.g. only one non-zero mass. For details see [17], [18]. To per- 



form the recursion numerically, it is important to use the multiple precision 
FORTRAN by D.Bailey (Jl9[) since tremendeous cancellations occur in this 
case. 

We have in some detail presented one approach for the calculation of 
the Taylor expansion of Feynman diagrams, others were worked out in Refs. 

@,|f]. The latter 

one is particularly suited for programming in terms of 
a formulae manipulating language like FORM. 



2.2. The method of analytic continuation 



Assume, the following Taylor expansion of a scalar diagram or a particu- 
lar amplitude is given C(pi,p2, . . .) = J2m=o a mV m = f(y) and the function 
on the r.h.s. has a cut for y > yQ. In the above case of H — > 77 one intro- 
duces y = q 2 /4m 2 with q 2 = {jp± — P2) 2 as adequate variable with yo = 1. 

The method of evaluation of the original series consists in a first step in 
a conformal mapping of the cut plane into the unit circle and secondly the 
reexpansion of the function under consideration into a power series w.r.t. 
the new conformal variable 



UJ = 1 -v i ~y/y\ (13) 

1 + Vi - y/yo 

By this conformal transformation, the y-plane, cut from yo to +00, is 
mapped into the unit circle (see FigJ^) and the cut itself is mapped on its 
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00 




Fig. 4. Conformal mapping of the y = g 2 /4rn^ -complex plane into the w-plane. 



boundary, the upper semicircle corresponding to the upper side of the cut. 
The origin goes into the point u = 0. 



After conformal transformation it is suggestive to improve the conver- 
gence of the new series w.r.t. u by applying the Pade method [ 22 ] , [ 23 1 . 
A convenient technique for the evaluation of Pade approximations is the 
e-algorithm of [22] which allows one to evaluate the Pade approximants 
recurrsively. 

Historically the first example considered was the two-loop three-point 
scalar ( planar) integral with the kinematics of the decay H — > 77: mj = 
and all other masses rrii = mt(i = 1,..,5). In this special case all Tay- 
lor coefficients can be expressed in terms of T -functions. Later a closed 
expression for arbitrary coefficients has been found in |2 



Table 1. Results on the cut (q 2 > 4m 2 ) in comparison with p4| . 





[14/14] 
Re 


Im 


Rcf.^ 
Re 


Im 


4.01 


11.935 


12.699 


11.9347(1) 


12.69675(8) 


4.10 


2.66245 


9.0955 


2.66246(2) 


9.0954(2) 


5.0 


- 1.985804823 


2.758626375 


- 1.98580(2) 


2.758625(2) 


10.0 


- 0.7569432708 


- 0.0615483234 


- 0.756943(1) 


- 0.061547(1) 


40.0 


- 0.045852780 


- 0.0645672604 


- 0.04585286(7) 


- 0.0645673(9) 


400.0 


+ 0.00008190 


- 0.0021670 


0.0000818974(3) 


- 0.002167005(3) 



Results for this kinematics on the cut are given in Table 1. The process 
H — > 77 was investigated before in Ref. |24j]. For the master integral under 
consideration in 1 24] all integrations but one could be performed analytically 
and only the last one had to be done numerically. Similarly, high precision 
is obtained on the cut in the approach of |ilf . 

Further examples of the efficiency of our method were given in Refs. [25[ 
and [26]. In the latter case a diagram of particular interest for the Z — » bb 
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decay has been considered, namely a planar diagram with a low threshold: 
mi = rri2 = = rriQ = m& = A.hGeV and 777.3 = m4 = = QIGeV . If 
this diagram is to be evaluated on the Z mass-shell, a precision of only some 
four decimals can be obtained with 30 Taylor coefficients, which is hardly 
good enough. On the other hand it is a very reasonable approximation in 
this case to put m& = from the very beginning and thus consider diagrams 
with zero thresholds, which should be evaluated more easily with higher 
precision. These will be considered in the next Sect. 



2.3. Two-loop vertex diagrams with zero thresholds 

Concerning the vertex diagrams, there are many different topologies 
contributing to a 3-point function in the SM. For our purpose of demonstra- 
ting the method, we confine ourselves to the planar case shown in Fig. ||a. 
Figs. |B|c,d presents infrared divergent diagrams. 



Pi P 2 



k 2 + Pl 

m 3 



k r k 2 




m 



1TL 



k 2 + P 2 
m 4 

k 1 + p 2 




a) generic b) c) d) 

Fig. 5. Planar diagrams with zero thresholds, a) — generic; b) — IR finite diagram; 
c),d) — IR divergent diagrams. 

As a typical (and very complicated case) we consider here Fig. ||b with 
equal non-zero masses |27|. There are two massless cuts so that we shall 
have the double logarithm in the expansion. The set of subgraphs in this 
case is given by 71 = V and the higher terms with 7 7^ V : 72 = {3456}, 73 = 
{1256}, and 74 = {56} (in curly brackets: the line numbers of the subgraphs 
7 in (^), line numbers are shown in Fig. ||). Note that 73 and 74 are dis- 
connected. After summing up all four contributions we see that the double 
and single poles in e cancel as well as the scale parameter fj,, with the result 
[x = q 2 /M 2 ) 



00 2 



n=0 j=0 



-XX 



{fo(x) + fi(x) H-x) + f 2 {x) ln 2 (-x)} , 



(14) 



12 



Table 2. Results for timelike q 2 for diagram 





[12/12] 




[15/15] 






Re 


Im 


Re 


Im 


0.05 


+2.948516245 


20.938528 


+2.948516245 


20.938528 


0.1 


-1.108116127 


16.04132127 


-1.108116127 


16.04132127 


0.5 


-4.820692281 


5.066080015 


-4.820692281 


5.066080015 


1.0 


-3.890154 


1.67549787 


-3.890156 


1.67549788 


1.5 


-2.904588 


0.42979 


-2.904581 


0.429778 


2.0 


-2.18294 


-0.06976 


-2.182981 


-0.069728 


10.0 


-0.191 


-0.208 


-0.194 


-0.215 



Table 3. The finite part for diagram 





[15/15] 




numerical results 




Re 


Im 


Re 


Im 


0.5 


81.17501719 


12.06458720 


81.1750 


12.0644 


1.0 


17.7659 


19.97799834 


17.7658 


19.9779 


2.0 


- 0.6047 


7.3759 


- 0.604 


7.376 


10.0 


- 0.572 


0.023 


- 0.570 


0.025 



where the fj n are now given in terms of rational numbers and C(2). f2( x ) 
can be summed analytically, yielding 

f 2 (x) = (In |1 + x\ - m0{l + x)) 2 /x 2 . (15) 

Thus, we have to Pade approximate fo and f\ only. Close to the second 
threshold at q 2 = M 2 the convergence is indeed excellent (see Table . 
It should be noted that for the physical application we have in mind, i.e. 
Z — > bb, this is just the case of interest. It is worthwhile to note the sharp 
increase for low q 2 due to ln 2 (— q 2 /M 2 ). 



Infrared divergent diagrams (Figs. ||c,d) have been successfully considered 
in 28]. For Fig. ||c the convergence is indeed excellent again. The different 
mass case, i.e. rriQ ^ mi = ?ti,2, had to be calculated numerically according 
to the procedure outlined in Sect. 2.1. For Fig. ||d the convergence is by far 
not that good: for q 2 /M 2 = 1 a precision of 6 decimals is still obtained with 
a [15/15] approximant, while for q 2 /M 2 ~ 10 only 2 decimals were obtained 
(see Table ||) where the "numerical results" were obtained by the method 
ofRef. 12T 



This much slower convergence implies that indeed it is very desirable to have 
at least in the case of only one non-zero mass a compact representation of 
the Taylor coefficients in terms of rational numbers, i.e. not to have to go 
through the mashinery of calculating the coefficients as described in Sect. 
2.1. Such a possibility will be described in the next section. 
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2-4- The Differential Equation Method 



We saw that to obtain the expansion of a diagram one has to go through 
a rather combersome machinery. The more coefficients are asked for, the 
more efforts and machine resourses are required. Thus it is very desirable 
to have analytic expressions for expansion coefficients whenever possible. 
This can be done with the aid of the Differential Equation Method (DEM) 
[j30f if only one non-zero mass occurs. The DEM allows one to get results 
for massive diagrams by reducing the problem to diagrams with simpler 
structure. 

Let us introduce a graphical notation for the scalar propagators (in 
euclidean space-time) 

1 _ a 1 



{7f ' tf + rr?) 

a and m refer to the index and mass of a line. Then one can derive the 

following recurrence relation for a massive triangle fjC| l 




(D — 2a\ — ct2 — 03) 



+«2 




-1 ai 




-2m 2 a\ 



-(ml+ml) / \ )+(a2^«3) 



along with some other graphical relations (see details in |30|| ). 

Using this technique we analysed in |}l]] the class of 3-point two-loop 
massive graphs. As an example for the diagram of Fig.|5|c we get 



v y 



(D - 4) J; 



5c 




4m z 




2m 2 ' 2 




(16) 



In the r.h.s. of (16) the last two terms can be combined, resulting in 
dJ/dm 2 while for the second term we proceed in turn as 



(D-A) 




(17) 



Thus we are left with simple diagrams (these can be done completely by 
Feynman parameters) and the derivative of the initial diagram w.r.t. m 2 . 
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The solution of the corresponding differential equation in terms of a series 
obtained from an integral representation reads 



T 2 (l+g) ~ (-x) n T 2 (n) 
5C = {q 2 f(m 2 t 2_, r(2n + l) 



i - ^(\n(x) + S l (n-l)j -^S 2 (n-1) 



1 5 1 

y<S?(n - 1) + 45i(n - l)5i(2n) - C(2) - ln(x)S 1 (n - 1) + - ln 2 (x) 



where 



1 

A similar formula was obtained for the diagram of Fig&I (see pi]]). 



2.5. Large mass expansion versus small momentum expansion 

If there are two or more different masses involved the coefficients a\ mn 
in (|j) are not just numbers any more but complicated functions of mass 
ratios. In this case one can try to perform a large mass (LM) expansion 
rather than a small momentum expansion. 

Here we just consider one example of a large mass expansion for the 
3-point function (a more detailed analysis will be given in fl3^1 ). Let us 
consider the diagram of Fig. ||a with two non-zero masses myy and mt op , 
contributing to Z decay in the SM. Thus we have q 2 = m 2 z . One can use the 
procedure described in Sect. 2.1. to expand this diagram in the external 
momentum squared. Then the actual parameter of expansion is q 2 /m 2 
and the radius of convergence is given by the lowest (non-zero) threshold 
i.e. q 2 = (mtop + m w) 2 ^ m %- Another possibility is to expand in the 
ratios q 2 /m, 2 op and rn^/m 2 ^, though one will see that this expansion does 
not work that well. 

The asymptotic expansion in the limit of the large mass is given by 
@. In our case there are 4 subgraphs that contribute, i.e. 71 = {123456}, 
72 = {126} 73 = {12345} and 74 = {12} (see Fig^a-d). By direct evaluation 
we find that there are induced poles of the order 1/e 3 in subgraphs while in 
the the sum they cancel which serves as a good check. For this particular 
diagram the result of the LM expansion looks like 

dia = £ A n = (l$ + 6«L + b^L 2 + 6(?) L 3 ) ( " (i 6) 



n=-l \ m top/ 
with L = log(m 2 op //i 2 ) and &r?'s being known functions of q 2 , m 2 ^ and fi 2 
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l- 5 --l i- i i 4 A- 

4 y crV *' :, v 

a) b) c) d) 

Fig. 6. A planar diagram a) with two different masses and subgraphs b),c) and 

d) contributing to the large mass expansion series. Dashed lines are massless; solid 
lines have light mass raw] bold lines have heavy mass m top . T stands for the Taylor 
expansion w.r.t. external momenta and light mass. 



In Table || we give the results of the numerical analysis for the graph at 
hand. We observed bad behavior of the series (|l6|): namely 10 coefficients 
give only 3 stable figures. Since the parameters of expansion q 2 /m^ op and 
both are of order 0.25 one would expect better convergence. We 



m 



w/ m top 



have found that this bad convergence remains true for some other graphs 
of the Zbb process [32|. Moreover no analytic properties are known in the 
1/m? variable and thus no conformal mapping like (13) is now available to 
improve the convergence. However we apply the Pade summation. Never- 
theless the LM expansion for 3-point functions in such a regime looks very 
attractive since it can be rather easily implemented in formulae manipula- 
ting languages like FORM even in the presence of more than one non-zero 
mass. In the worst case one achives 0.1% accuracy with 10 coefflciients. 
Our programs allow us to get 15-20 coefficients in a reasonable time. 

Table 4. Terms in the LM expansion corresponding to formula (|l^) at [i = mtop = 
180 , m\v = 80 , q 2 — 90 2 . The last two lines give Pade approximants obtained with 
10 terms of the LM expansion and 8 Taylor coefficients of the small q 2 expansion. 



n 


Re An 


Im A„ 


-1 


-6.35040 


+ 29.98705 





-3.85690 


- 6.19256 


1 


-1.74681 


- 4.34057 


2 


0.22329 


- 1.30291 


3 


0.65892 


- 0.21984 


4 


0.47058 


+ 0.00121 


5 


0.25376 


+ 0.01435 


6 


0.13099 


+ 0.00492 


7 


0.07311 


+ 0.00073 


8 


0.04517 


- 0.00009 


9 


0.02984 


- 0.00009 


10 


0.02042 


- 0.00002 


LM 


-9.996 


17.9527 


small-g 


-9.9926682590236 


17.95215366130182 



2.6. Conclusion 
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Any involved calculation in field theory nessesarily consists of two parts: 
l)automatic generation of Feynman diagrams and source codes and 2)tech- 
niques of evaluating scalar Feynman diagrams. In this paper we presented 
both. Still some efforts have to be applied to handle the numerators of 
diagrams. We are working also on the automation of the renormalization 
procedure. 
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